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ALIAS-FREE  WIGNER  DISTRIBUTION  FUNCTION  AND 
COMPLEX  AMBIGUITY  FUNCTION  FOR  DISCRETE-TIME  SAMPLES 

INTRODUCTION 


The  attributes  of  the  Wlgner  distribution  function  (WDF)  have  come  under 
close  scrutiny  In  recent  years;  see,  for  example,  [1,2,3]  and  the  references 
listed  therein.  However,  the  numerical  calculation  of  the  WDF  from  discrete 
time  data  still  suffers  from  the  belief  that  the  sampling  rate  of  a  given 
time  waveform  must  be  twice  as  large  for  computation  of  an  alias-free  WDF, 
as  the  rate  required  for  reconstruction  of  the  original  continuous 
waveform.  If  true,  this  would  double  the  number  of  data  points  that  must  be 
collected  to  cover  a  given  time  interval,  and  greatly  increase  the  number  of 
subsequent  computations.  This  contention  applies  to  the  complex  analytic 
waveform  as  well  as  to  a  real  waveform. 

It  Is  the  purpose  of  this  report  to  establish  the  fact  that  the  sampling 
rate  need  not  be  doubled,  and  that  an  alias-free  WDF,  as  well  as  complex 
ambiguity  function  (CAF),  can  still  be  quickly  and  efficiently  obtained, 
provided  that  fill  the  Information  In  the  available  data  stream  Is  extracted 
and  properly  processed.  Some  recent  effort  on  this  topic  [4,5,6]  did  not 
discover  the  particular  complete  set  of  processing  required,  leading  to  the 
conjecture  [5,  page  1068]  that  It  was  not  possible  to  accomplish  the  desired 
goal  for  the  WDF. 
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We  will  show  not  only  that  the  desired  goal  can  be  achieved,  but  that 
the  required  data  processing  for  an  alias-free  WDF  and  CAF  Is  strikingly 
simple.  Our  approach  to  the  solution  Initially  Involves  the  four 
time/frequency  domains  associated  with  the  WDF  and  Its  various'  Fourier 
transforms.  However,  in  hindsight,  an  extremely  simple  and  direct  method  of 
obtaining  the  WDF  and  CAF  will  be  presented,  which  requires  only  FFTs  (fast 
Fourier  transforms)  for  Its  Implementation. 

An  alias-free  discrete  WDF  and  CAF  have  been  achieved  In  [7]  and  [8],  by 
means  of  Interpolating  either  the  waveform  time  samples  or  the  spectrum 
frequency  samples.  Also,  the  ranges  and  required  Increment  sizes  In  time 
and  frequency  of  the  various  two-dimensional  functions  have  been  carefully 
scrutinized  In  [7],  by  discrete  Fourier  transform  techniques.  However,  that 
approach  does  not  Illuminate  how  the  various  aliasing  lobes  Interact  and  can 
be  controlled.  Furthermore,  we  utilize  a  continuous  approximation  approach 
(rather  than  a  discrete  Fourier  technique),  which  lends  tremendous  Insight 
Into  the  shortcomings  of  current  processing  methods  and  brings  out  the 
fundamental  properties  of  the  various  two-dimensional  functions  and  their 
domains  of  definition.  The  final  discretization  In  time  and  frequency  Is 
only  done  with  deference  to  practical  computer  evaluation. 

In  fact,  we  will  not  define  a  discrete  WDF  or  CAF  here.  Instead,  we 
attempt  to  recover  the  WDF  and  CAF  of  the  original  continuous  time  waveform, 
by  developing  approximations  and  then  controlling  or  eliminating  the  errors 
In  these  approximations.  Only  after  this  is  accomplished,  do  we  then  address 
discretization  of  the  time  and  frequency  arguments  of  the  two-dimensional 
functions  of  Interest. 
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NOTATION 


For  economy  of  presentation,  a  number  of  notatlonal  and  manipulative 
shortcuts  are  employed  here.  We  have  collected  them  all  together  at  this 
point,  and  will  employ  them  freely  later,  with  minimal  comment.  We  define 


rect(x) 


1  for  |x\  <  1/2 
0  otherwise 


si nc( x)  -  )  for  all  x  . 


(1) 


(2) 


The  symbols  j*  and  without  limits  denote  that  Integration  and  summation 
are  to  be  conducted  over  the  complete  range  of  nonzero  Integrand  and  summand, 
respectively. 


The  convolution  of  two  functions  g(x)  and  h(x)  is  denoted  by 

g(x)«  h(x)  -  Jdu  g(u)  h(x  -  u)  .  (3) 

The  two-dimensional  convolution  of  two  functions  Is 

xy  rr 

g(x,y)  ®  h(x,y)  -  JJ  du  dv  g(u,v)  h(x  -  u,y  -  v)  .  (4) 

The  Fourier  transform  of  a  time  domain  function  s(t)  Into  Its  spectrum 
In  the  frequency  domain  f  Is  according  to  the  pair  of  relations 

S(f)  ■  j*  dt  exp(-12irft)  s(t)  , 

s ( t )  -  pdf  exp(  1 2irf t)  S(f)  .  (5) 
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Then  the  Fourier  transform  of  a  product  of  time  functions  is  equal  to  the 
convolution  of  the  two  spectra: 

Jdt  exp(-i2irft)  a(t)  b(t)  =  A(f)$B(f>  .  (6) 

The  infinite  impulse  train  in  time  t,  with  spacing  a,  is 

« 

a  ^4(t  -  na)  ,  n  Integer  .  (7) 

n 

Its  Fourier  transform  is  another  infinite  Impulse  train  In  frequency  f,  with 
reciprocal  spacing: 

j*  dt  exp(-12irft)  a  4(t  -  na)  -  ^a(f  -  J)  .  (8) 

n  n 

Combination  of  (6)  and  (8)  leads  to  a  very  useful  relation  that  Is  employed 
frequently  In  the  following: 

j'  dt  exp( -1 2irf t )  a(t)  a  ^4(t  -  na)  » 

n 

-  A(f)  e  -  J)  -  2  -  J)  •  O) 

n  n 

The  discrete  Fourier  transform  operation  arises  frequently;  consider 

Z(n)  5  exp(-12irkn/N)  z(k)  for  all  n  .  (10) 

k 
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Ihe  periodicity  of  Z(n)  means  that  it  only  need  be  computed  for  one  period, 
namely  0  <  n  <  N  -  1.  The  absence  of  limits  on  the  sum  In  (10)  means  that 
it  goes  from  k  -  -®  to  +<».  However,  since  z(k),  z(k  t  N),  z(k  ±  2N), 

...  all  receive  the  same  weight  in  (10),  regardless  of  the  value  of  n,  the 
values  of  {z(kjj  can  be  "collapsed1'  according  to 


z(k) 


z(k  ♦  .IN) 


for  0  i  k  <  N 


otherwise 


1 

► 


01) 


and  (10)  becomes  Identically 


N-l 

Z(n)  «  ^  exp(-12irkn/N)  z(k)  for  all  n  .  (12) 

k-0 

For  n  highly  composite,  FFT  routines  can  be  employed  for  efficient  evaluation 
of  (12)  for  0  <  n  <  N  -  1,  The  manipulation  of  (10)  Into  (12)  Is  called 
collapsing  (or  prealiasing),  and  the  operation  In  (11)  is  modulo  N  addition. 
The  nonzero  values  of  {z^  In  (10)  can  occur  anywhere  on  the  k-axls,  and 
there  can  be  an  arbitrary  number  of  them;  nevertheless,  (12)  is  an  Identity 
with  (10).  The  value  of  Z(n)  for  any  n  can  be  obtained  Immediately  from  the 
FFT  output,  by  looking  up  the  value  in  location  n  modulo  N. 
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WAVEFORM  CHARACTERUAT ION  AND  RATE  OF  VARIATION 

The  continuous  complex  waveform  of  Interest  Is  s(t),  with 
F*uri*r  spectrum 

$(f)  «•  J  dt  exp{ -12irft)  s ( t )  for  all  f  .  (13) 

We  presume  here  that  spectrum  S(f)  Is  band  11ml ted (  with  total  extent  F  Hz; 

1  .e. , 

S(f)  -  0  for  \f\  >  F/2  .  (14) 

Notice  In  figure  1  that  spectrum  S(f)  Is  centered  at  f  -  0,  without  loss  of 
generality,  since  waveform  s(t)  could  be  multiplied  by  exp(-12irfQt)  to 
downshift  It  by  f  Hz,  to  any  convenient  center  frequency,  as  desired. 

If  we  were  given  a  real  waveform,  we  would  replace  It  by  Its  analytic 
waveform  or  complex  envelope,  thereby  allowing  the  minimal  possible 
time-sampling  rate  that  can  still  exactly  represent  and  recover  the  complex 
waveform.  This  sampling  rate  is  half  that  required  for  sampling  the 
corresponding  real  lowpass  waveform,  without  loss  of  Information. 
Nevertheless,  the  decreased  sampling  rate  applied  to  the  complex  waveform  Is 
still  sufficient  to  get  an  alias-free  WOF  and  CAF.  (Of  course,  the  samples 
are  now  complex,  whereas  they  were  formerly  real  for  the  real  waveform  case.) 
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SPEC  I  RUM  OF  SAMPLED  WAVEFORM 


Waveform  s(t)  Is  sampled  at  time  Increment  4  seconds,  yielding  samples 

{s(ka)^  for  all  Integer  k  .  (15) 

The  spectrum  of  this  sampled  waveform  Is  defined  by  means  of  a  Trapezoidal 
approximation  to  defining  Integral  (13): 


5(f)  =  A  ^ 


exp(-12*fAk)  s(kA)  * 


-  f  dt  exp( -12irf t)  s(t)  A  «(t  -  kA)  - 
J  k 

S(f)  $  2A(f  “  A}  “  ZS(f  “  A>  f°r  311  f  ' 


(16) 


where  we  used  (6)-(9).  The  approximating  spectrum  3(f)  has  period  l/A  In  f 
and  Is  depicted  In  figure  2.  It  will  have  nonoverlapping  aliasing  lobes  If 


(17) 


This  fundamental  sampling  rate  condition  will  be  presumed  to  be  true, 
henceforth.  In  fact,  In  order  to  keep  the  number  of  samples  {s(kA)]  small, 
(17)  will  be  presumed  to  be  closely  met.  It  Is  very  Important  to  minimize 
the  number  of  samples  that  must  be  manipulated,  so  that  the  computational 
burden  In  evaluating  the  WDF  Is  not  overwhelming. 


Another  Interpretation  of  approximation  S(f)  Is  afforded  by  line  2  of 
(16):  ?(f)  Is  the  spectrum  (Fourier  transform)  of  the  signal  s(t)  sampled 
(multiplied)  by  the  Infinite  Impulse  train  at  spacing  a.  This  alternative 
Interpretation  will  also  arise  later,  when  we  Investigate  sampling  relative 
to  the  WDF  and  Its  various  two-dimensional  transform  domains. 
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Since  sampling  rate  condition  (17)  is  presumed  to  be  met,  then  spectrum 
S(f )  In  figure  2  can  be  gated  with  a  rectangular  function,  and  S(f)  can  be 
recovered;  l.e., 

S(f)  «  S(f)  rect(fa)  for  .all  f  .  (18) 

Therefore,  the  time  waveform  s(t)  can  also  be  recovered  exactly,  for  all  t, 
by  means  of  Inverse  transform  (5); 

s(t)  ■  j*  df  exp(12irft)  S(f)  rect(fA)  - 

“  s(t)  %  \  slnc(J)  .  (19) 

However,  since  from  line  2  of  (16),  product  waveform 

s(t)  -  s(t)  &  "  kA)  “  A  2*(kA)  -  kA)  *  (20) 

k  k 

then  (19)  becomes 

s(t)  -  2s(kd)  s1nc(^  -  k)  for  all  t  ,  (21) 

k 

which  Is  the  standard  Interpolation  formula  for  a  bandllmited  waveform. 

It  should  be  pointed  out  here  that  (21)  is  not  an  attractive 
computational  procedure,  and  that  an  excellent  alternative  Is  available. 
Namely,  from  (16),  compute  from  the  available  samples, 

S(f)  -  A  exp(-12*fAk)  s(kA)  for|f|<|J,  (22) 

k 

and  then  use  the  top  line  of  (19)  to  recover  waveform 
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_1 

2a 

s(t)  -  df  exp(l2irft)  S(f)  for  all  t  .  (23) 

zl 

2a 


The  reason  this  gating  procedure  Is  attractive  Is  that  (22)  and  (23)  can  both 
be  done  by  FFT  procedures.*  Also,  this  seemingly  trivial  sidelight  will 
reoccur  in  WDF  and  CAF  reconstruction,  where  It  will  have  a  significant 
Impact. 


Notice  that  we  have  not  defined  a  discrete  spectrum,  per  se.  Rather,  we 
have  concentrated  on  getting  an  approximation  ?(f)  to  the  original  continuous 
spectrum  S(f),  both  defined  for  all  f.  If  sampling  condition  (17)  Is  met, 
a  <  1/F,  the  approximation  affords  the  possibility  of  exact  recovery  of  S(f) 
at  any  f.  This  philosophy,  namely  avoiding  arbitrary  definitions  of  discrete 
functions,  In  favor  of  direct  approximations  to  the  desired  continuous 
functions,  Is  pursued  throughout  this  report.  It  Is  believed  that  this 
clarifies  the  fundamental  limitations  and  processing  that  must  be  performed 
In  order  to  achieve  the  desired  quantities.  Finally,  after  demonstrating 
the  viability  of  this  approach,  In  order  to  reduce  the  mathematical  equations 
to  practical  calculations,  we  discretize  the  time  and/or  frequency  arguments 
of  the  approximations,  as  appropriate,  and  manipulate  the  equations  Into 
attractive  FFT  forms.  We  end  up,  of  course,  with  discrete  data  processing 
forms  that  are  suitable  for  efficient  computer  realization,  but  the 


*Actually,  termination  of  the  sum  In  (22)  at  finite  k  limits  will  yield  an 
approximation  to  5(f);  the  error  can  be  controlled  to  any  desired  degree  by 
taking  enough  terms.  Also,  the  Integral  In  (23)  will  have  to  be 
approximated,  say,  by  the  Trapezoidal  rule;  the  attendant  time-aliasing  can 
be  minimized  by  choosing  the  frequency  Increment  small  enough.  These  details 
will  be  Investigated  later. 
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discretization  in  time/frequency  Is  deferred  to  the  latest  possible 
location,  since  it  is  not  fundamental  to  the  ideas  of  controlling  or 
eliminating  aliasing. 

The  sampling  Increment  a  will  not  be  set  equal  to  1  in  this  report,  for 
several  reasons.  It  Is  easier  to  keep  track  of  dimensions,  and  dimensional 
checks  on  the  equations  are  accomplished  more  readily.  It  is  also  easier  to 
obtain  physical  Interpretation  of  time  instants  and  Increments,  as  well  as 
frequency  limits  and  bandwldths.  Finally,  it  will  be  seen  to  eliminate 
confusion  and  ambiguity  as  to  precisely  where  time  and  frequency  samples  of 
the  temporal  correlation  function,  WDF,  and  CAF  are  being  taken;  the 
Importance  of  this  last  point  can  not  be  overemphasized. 

EXAMPLES 


It  Is  very  Informative  at  this  point  to  consider  a  couple  of  continuous 
waveforms  and  their  corresponding  WDFs,  In  terms  of  their  rates  of  variation. 
Consider  first,  spectrum 


S(f) 


(24) 


for  which  the  waveform  Is 


s(t)  «•  slnc(Ft)  -  for  all  t  . 

The  corresponding  WOF,  at  time  t  and  frequency  f,  Is  [9,  (10)] 


(25) 
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W(t.f)  r  J  dT  exp(-i2irft)  s ( t  +  |)  s* (t  -  |)  * 

■  j*  dv  exp( I2trwt)  S(f  +  |)  S*(f  -  |)  - 

for  |f)  <  |  j 

>  for  all  t 

*  '  2  -^“)s1nc  C2Ft(l  “  2^)]  rect(£)  ,  (26) 

Then,  for  Instance,  the  slice  of  the  WOP  at  zero  frequency, 

W(t,0)  -  jf  s Inc ( 2Ft )  ,  (27) 

varies  as  fast  as  waveform  $(t)  In  (25).  Therefore,  although  sampling 
s(t)  In  (25)  with  time  Increment  A  <  1/F  Is  sufficient  to  reconstruct  s(t), 
we  need  a  time  Increment  half  as  large  In  order  to  adequately  sample  slice 
(27)  of  the  WDF  at  f  -  0.  In  fact,  W(t,f)  In  (26)  varies  faster  with  t  th«n 
s(t)  does,  whenever  if  1  <  F/4. 

This  example  points  out  that  the  WDF  must  be  computed  twice  as  finely  as 
the  waveform  samples,  lest  Important  Information  about  the  energy 
distribution  of  s(t)  in  t,f  space  be  lost.  In  fact,  If  (27)  were  computed 
at  time  points 

tn  -  (n  +  ^)  £  for  all  n  ,  (28) 
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which  have  time  increment  A  -  1/F,  then  the  WDF  values  obtained  would  be 

w((n  f  J)  J-  ,  o)  -  jj-  sinc(2n  ♦  1)  •  0  for  all  n  .  (29) 

We  would  be  led  to  believe  from  samples  (29)  that  there  Is  no  energy  along 
the  f  *  0  line  In  t,f  space,  whereas  continuous  version  (27)  Indicates  a 
considerable  contribution. 


(30) 


(31) 


This  WOP  varies  faster  with  t  than  s(t)  does,  and  faster  with  f  than  S(f) 
does.  In  fact,  the  rates  of  variation  of  W(t,f)  and  |i(t)|  2  are  the  same 
with  t,  while  those  of  W(t,f)  and  Js(f)|?  are  the  same  In  f. 


Both  of  the  examples  above  Illustrate  the  need  to  compute  the  WDF  at 
finer  Increments  than  are  adequate  for  the  time  waveform  or  spectrum. 
However,  this  does  oat  mean  that  the  time  waveform  need  be  sampled  more 
frequently  than  requirement  (17).  Rate  (17)  Is  fine  for  sampling  waveform 
s(t),  but  the  corresponding  WDF  can  and  must  then  be  computed  at  finer 
increments. 
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TWO-DIMENSIONAL  CONTINUOUS  FUNCTIONS 

For  a  continuous  waveform  s < t )  with  spectrum  S(f),  there  are  four  useful 
two-dimensional  characterizations.  The  first  Is  the  continuous  temporal 
correlation  function  (TCF) 

R ( t , t )  -  s(t  +  J)  s*(t  -  |)  for  all  t,t  .  (32) 

Variable  t  Is  absolute  time  In  seconds,  while  t  is  relative  time  or  time 
separation.  The  corresponding  Wlgner  distribution  function  (WDF)  Is  a 
Fourier  transform  on  t: 

W(t,f)  -  j*  dr  exp( -12irf t)  R(t,t)  for  all  t,f  .  (33) 

The  alternative  Fourier  transform  on  t  yields  the  complex  ambiguity  function 
(CAF): 

x(v,t)  ■  J  dt  exp(-12tvt)  R(t,T)  for  all  v,t  .  (34) 

Functions  W  and  x  are  two-dimensional  Fourier  transforms  of  each  other. 

Finally,  completing  both  routes  (by  t  or  by  0,  we  have  the  spectral 
correlation  function  (SCF)  as  another  Fourier  transform,  according  to 
several  equivalent  forms 
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A(v.f) 


J*  dt  exp(  -I2irut)  W(t,f )  - 
J dt  exp( -1 2irf t)  x(v,t)  ■ 


■  dt  dt  exp(-12irvt  -  1 2irf t )  R(t,t)  ■ 

-  S(f  +  j)  sV  -  j)  for  all  v,f  .  (35) 

This  last  relation  In  terms  of  the  spectrum  S(f)  of  waveform  s(t)  will  turn 
out  to  be  extremely  Important  and  useful.  It  also  enables  Interpretation  of 
f  as  absolute  frequency  In  Hz,  while  v  Is  relative  frequency  or  frequency 
separation. 

The  names  for  the  TCP  and  SCF  have  been  drawn  from  the  similarity  of 
their  forms  In  (32)  and  (35),  respectively,  to  correlation  operations.  The 
latter  name  Is  also  used  In  (10,  (5)-(7)]  for  a  similar  quantity. 

Recalling  the  bandllmlted  character  of  S(f)  In  (14)  and  figure  1,  we  see 
that  SCF  A(v,f)  In  (35)  can  be  nonzero  only  when 

j f  *  2)  <  f  •  (36) 

This  region  In  the  two-frequency  domain  (M  plane)  Is  depicted  In  figure 
3.  It  Is  a  diamond-shaped  region  centered  at  the  origin  of  the  v,f  plane. 
Outside  this  diamond,  SCF  A( v , f )  Is  'identically  zero.  Thus,  a  bandllmlted 
spectrum  S(f)  Is  reflected  In  the  «,f  domain  as  a  diamond-limited  SCF  A(v,f). 
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Since  (35)  can  be  inverted  to  give 

W(t,f)  *  j* dv  exp(l2irvt)  A(v,f)  , 

X<v,t)  -  [  df  exp(i2irft)  A(v,f )  ,  (37) 

It  follows  from  figure  3  that 

W(t,f)  -  0  for  |f|  >  F/2  , 

X(v,r)  *  0  for  |v|  >  F  .  (38) 

Thus  the  WDF  and  CAF  are  band  1 1  ml  ted  In  their  respective  frequency  variables. 
These  properties  will  be  useful  later  when  we  study  the  effects  of  aliasing 
In  the  various  domains.  More  generally,  the  extents  and  rates  of  variation 
of  the  TCF,  WDF,  CAF,  and  SCF  are  summarized  In  appendix  A. 

The  following  symmetry  properties  on  the  TCF  and  SCF  reduce  computational 
effort  by  a  factor  of  two! 

R(t,~r)  -  R*(t,t)  , 

A(-v,f)  -  Aw(«,f )  . 

These  follow  Immediately  from  (32)  and  (35),  respectively. 
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TWO-DIMENSIONAL  FUNCTIONS  FOR  DISCRETE -TIME  SAMPLES 


The  available  data  samples  of  waveform  s(t)  are,  as  given  1i  i15), 


js(kA)}  for  Integer  k  . 
SAMPLED  TEMPORAL  CORRELATION  FUNCTION 


From  these  values,  the  totality  of  information,  that  can  be  computed 
regarding  the  continuous  TCF  R(t , t)  In  (32),  are  the  sets  of  discrete 
values 


R(mA,2qA)  -  s(mA  +  qA)  s  (mA  -  q A) 


and 


^m  +  J)A,  (2q  +  1)a)  -  s(mA  +  qA  +  A)  s*(mA  - 


^for  Integers 
m  and  q, 


qA), 


(39a) 


(39b) 


Thus  both  the  t  and  t  variables  In  R(t,?)  are  discretized,  as  Indicated 
In  figure  4.  However,  observe  that  the  available  Information  is 
interspersed  In  the  t,t  plane.  Thus,  for  fixed  t,  the  separation  In 
available  t  values  Is  2A,  not  A,  whether  t/A  Is  Integer  or  half-integer. 
Similarly,  for  fixed  t,  the  separation  In  available  t  values  Is  a,  not 
A/2,  whether  t/A  Is  an  even  or  odd  Integer.  This  lack  of  intermediate 
values  in  both  slices  1$  what  has  led  (in  the  past)  to  Incomplete  processing 
of  the  available  information.  What  is  needed  to  solve  the  aliasing  problem 
Is  a  combination  of  ill  the  Interspersed  Information  In  figure  4  Into  a 
single  unified  two-dimensional  description.  That  solution  will  be  found  to 
reside  in  the  SCF  domain,  v,f. 
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Figure  4.  Available  Values  of  TCF  R(t,t) 


APPROXIMATE  WIQNER  DISTRIBUTION  FUNCTION 

Guided  by  definition  (33),  we  adopt  the  following  Trapezoidal 
approximation  to  the  continuous  WOF  W(t,f)  at  time  t  «  m,  m  Integer: 

Wfl(ma,f)  =  2a^>  exp(-l2irf2qa)  R(ma,2qa)  for  all  f  .  (40) 

d 

Subscript  a  on  Wfl  denotes  that  It  Is  only  an  approximation  to  the  true 
continuous  W.  Notice  that  the  t  Increment  In  (40)  Is  2A,  as  It  must  be, 
according  to  (39a)  and  figure  4;  we  are  taking  a  vertical  slice  at  t  »  ma  In 
figure  4.  The  approximation  In  (40)  Is  always  real.  It  utilizes  only  the 
upper  line  of  Information  available  In  (39).  Notice  also  that  this  function 
Is  defined  for  all  f. 

However,  there  Is  an  additional  approximation  to  W(t,f)  available  at 
time  t  «  (m  +■  ^)a  ,  by  use  of  the  bottom  line  of  (39);  namely,  guided  again 
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by  definition  (33),  we  have  Trapezoidal  approximation 


Wa  +  2^,f)  s  2A  ^  exp[-i2irf (2q  +  1)A]  R^(m  +  ^)A,  (2q  +  1)A^  for  all  f  . 


m  is  still  an  integer.  Again,  the  t  Increment  is  2A,  but  it  is  shifted  by 
A,  as  (39b)  and  figure  4  dictate.  We  are  now  taking  a  vertical  slice  at 

t  ■  (m  +  ^)a  In  figure  4;  this  is  in  keeping  with  the  philosophy  developed 

earlier  in  ( 24) -*( 31 ) .  Approximation  (41)  is  also  real. 


Equations  (40)  and  (41)  can  be  developed  into  some  Informative  forms; 
from  (40)  and  (9),  approximation 

Wa(mA,f)  -  J dt  exp(-12irfT)  R(mA,t)  2A^T  4(t  -  2qA)  *= 

q 

-  W(mA,f)  9  4(f  -  2fi)  ‘ 

q 

-  W(mA, f  -  2^)  for  all  f  .  (42) 

q 

Similarly,  (41 )  yields 


Wa(jm  *  2)A,f)  "  J  dT  exP(-12fffT)  R($m  +  2)a,T)  U  J 

q 

-  w(jm  ♦  £)A,f)«  ^>(-l)q  W  -  m 

q 

*^>(-1)^  W  ^m  +  |^)A,  f  -  2^  for  all  f  . 


-  2qA  -  A)  * 
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The  (-1)q  factor  is  due  to  the  time  delay  of  a,  seconds  in  the  impulse 
train;  more  generally,  for  delay  tq, 

exp(  -1 2irf t )  2A^jT  4(t  ~  2qA  -  tq)  ■ 

q 

-  exp(-12«fr0>  ^  4{f  -  2?)  * 

q 

*^>exp(-1irq  tq/A)  4(f  -  .  (44) 

q 

The  two  relations,  (42)  and  (43),  are  equivalent  to  those  given  in 
[4,  (9)]  and  [6,  (14)].  Observe  that  the  aliasing  lobes  are  separated  by 
only  (2a)  1  in  the  frequency  domain,  not  1/A  as  was  the  case  for  the  signal 
aliasing  lobes  In  figure  2.  These  approximations,  Wa,  are  Illustrated  In 
figure  5,  for  two  adjacent  time  Instants  at  mA  and  (m  *  |)A  .  We  have  used 
property  (38)  In  drawing  figure  5. 

In  order  for  either  approximation,  by  Itself,  to  be  free  of  aliasing,  we 
would  need 

I  <  2A  "  2  '  1,e**  4  <  2?  * 

This  relation,  obtained  directly  from  both  plots  In  figure  5,  Is  the  usual 

it 

one  quoted  regarding  an  alias-free  WOF.  It  Is  seen  to  require  a  sampling 
rate  twice  as  fine  as  (17).  If  we  satisfy  (17),  but  not  (45),  then  the 
approximations  in  figure  5  are  significantly  aliased. 

The  case  where  s(t)  Is  a  real  waveform  Is  treated  In  appendix  A. 
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The  two  approximations,  (40)  and  (41),  use  all  the  available  information 
(39)  about  the  TCF  R(t,t)  In  the  t.t  plane.  However,  we  cannot  average 
these  two  WDF  approximations,  In  hopes  of  cancelling  out  the  close-in  lobes 

centered  at  f  *-  t(2A)-1,  because  the  two  times,  mA  and  (m  f  g)A,  are  not 
identical.  (This  timing  observation  Is  one  reason  for  keeping  A  Itself  In 
all  the  equations,  rather  than  setting  a  equal  to  1  and  losing  track  of  the 

meaning  of  m  vs.  m  +  g) •  Nor  can  we  discard  either  one  of  approximations 

(42)  and  (43),  especially  if  criterion  A  <  1/F  is  closely  met;  the  examples 
in  (24) -(31)  amply  demonstrate  the  rapid  variation  of  the  WDF  with  time  t. 

We  see  from  (40)  and  (41)  that  approximations  to  the  continuous  WDF  are 
available  at  discrete  time  values  with  separation  A^  *■  A/2  and  at  a 
continuum  of  f  values.  Thus  we  have  succeeded  In  eliminating  the  discrete 
nature  of  one  of  the  two  initial  time  variables  in  the  1CF,  namely  t. 
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APPROXIMATE  SPEC  I KAL  CORRFLATIQN  FUNCTION  VIA  WUF 

Guided  by  line  1  of  definition  (35),  we  adopt  the  following  Trapezoidal 
approximation  to  the  original  SCF  A(v,f): 

A<w)(v,f)  2  f  2.  ®xp(i2™  *2)  wa(fi2»f)  “ 
n 

^  +■  j  exp(-12irW  Wa(Q|,f)  for  all  v,f  .  (46) 

In  even  n  odd  J 

Notice  that  this  function  Is  defined  on  a  continuum  in  v,f  space.  The 
superscript  w  on  approximation  A  denotes  the  fact  that  we  have  used  the 

Q 

WOF  route  to  get  into  the  v,f  plane.  The  Increment  on  t  In  (46)  Is  At  «  a/2, 
In  keeping  with  the  available  WDF  values  In  (40)  and  (41)  together.  This 

fwl  (wl  * 

approximation  satisfies  the  symmetry  rule,  Aj  '(-v,f)  ■  A^  '(v,f)  ,  just  as 

for  the  original  SCF  A(u,f). 

We  let  n  *  2m  in  the  first  sum  in  (46),  and  let  n  ■  2m  +  1  In  the  second 
sum.  There  follows,  upon  use  of  (42)  and  (43), 

A^w^(v,f)  ■  g  exp(  -12ir«am)  Wfl(ma,f)  + 
m 

+  ^exp^-12irva(m  +  ^W^m  +  “ 

m 

*  2  ]>  exp(-12irvam)  W(ma,f  -  +• 

q  m 

+  exp[-12i»va(m  +  ^)]  W^m  f  £)a,f  -  ^  •  (47) 

q  m 
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while 


A  ^  exp(  -12ir«Am)  W(mA,f')  » 
m 

j*  dt  exp(  -12irvt)  W(t.f')  A  6(t  -  m  A)  - 

in 

-  A(w,f')e  ^  *(v  -  J)  - 

m 

-  ;,f ) , 


a]>  exp[-l2irvA(n>  +  ^)]  W +  |)A,f^ 


J  dt  exp(- 


12irut)  W(t.f')  A 


24- 


A(v.f')  ®2(“1)m  i(w  “  A) 


^(-Dm  A( v  -  ®,f ') 
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The  use  of  these  two  relations  In  (47)  yields  the  approximate  SCF 


a‘w><».o  ■  1 2  2 A("  -  ?• f  -  - 

(50a) 

q  m 

»<•  -?• f  -j?)  * 

(50b) 

q  m 

2  2  A<»  -  !•  f  -  »})  for  a11  v»f  • 

(51) 

q  m 

q-Hti  even 

We  now  have  a  function  defined  on  a  two-dimensional  continuum  In  the 
two-frequency  domain,  v,f. 

At  this  point,  the  reason  for  pursuing  the  use  of  all  the  available 
Information  becomes  obvious.  All  the  close-in  lobes  that  caused  problems 
have  precisely  cancelled  In  the  SCF  domain!  Figure  6  depicts  the  regions  In 
the  v,f  plane  where  approximation  A^fv.f)  In  (51)  can  be  nonzero;  see 

Q 

also  figure  3. 

The  SCF  term  In  (50a)  corresponds  to  use  of  only  the  Information  about 
the  TCF  given  In  (39a),  within  a  factor  of  2,  whereas  (50b)  arises  from 
(39b).  Term  (50a)  by  Itself  contains  all  the  aliasing  lobes  centered  at 

v  -  ®,f  ■  2^,  with  separations  Av  -  1/a,  af  »  (2a)"1;  condition  (17)  Is  then 
Insufficient  to  prevent  overlap,  and  (50a)  Is  seriously  aliased.  A  similar 
situation  exists  for  (50b)  by  Itself.  It  is  only  the  average  of  these  two 
pieces  of  Information  that  succeeds  In  elimination  of  the  troublesome 
close-in  aliasing  lobes  In  v,f  space.  And  It  Is  only  in  this  last  domain, 
where  the  functions  are  continuous  In  both  variables,  that  this  average  can 
be  conducted. 
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There  will  be  no  overlap  of  any  of  the  remaining  aliasing  lobes  In 
figure  6  If  we  choose,  as  In  (17), 

A  <  £  .  (52) 

Notice  that  we  do  have  to  require  a  <  (2F)-1,  In  order  to  avoid  the 
overlap.  Furthermore,  if  we  define  the  diamond  gating  function  (see 
figure  6) 


then  we  can  recover  exactly  the  original  SCF  from  approximation  (51) 
according  to 

A<w)(v,f)  D(v,f )  -  A(v,f )  for  all  «,f  ,  (54) 

tut  recovery  of  A(v,f)  Is  tantamount  to  recovering  the  exact  continuous  WDF, 
since 


W(t,f)  -  jdv  exp(12irvt)  A(v,f )  for  all  t,f  .  (55) 

Thus,  criterion  (52)  Is  sufficient  to  guarantee  the  possibility  of  getting 
an  alias-free  WDF  from  discrete-time  data. 

Additional  Interpretations  of  (54)  and  a  simple  method  of  computing 
approximation  A^fu.f)  are  addressed  in  the  next  section,  after  we 
have  also  looked  at  the  route  to  the  v,f  plane  by  way  of  the  CAF. 
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APPROX IMA1 t  CUMPLtX  AMBIGUITY  FUNCIION 

Based  on  definition  (34),  we  utilize  the  Trapezoidal  approximation  to 
the  continuous  CAF  x(u,t)  at  delay  t  «=  2qA,  q  Integer: 

Xa(w,2qa)  =  A^  exp(~12irvAm)  R(mA,2qA)  « 
m 

■  Jdt  exp( -I2irvt)  R(t,2qA)  A  4(t  -  mA)  - 

m 

■  x(«,2qA)  4(v  -  a)  ■ 

m 

■  ^>  x(*  -  S*  2£1A)  for  a11  w  •  (56) 

m 

Notice  that  the  t  Increment  Is  A,  as  It  must  be,  according  to  (39a)  and 
figure  4;  we  are  taking  a  horizontal  slice  at  t  «  2qA  In  figure  4. 


However,  there  Is  an  additional  approximation  available  to  the  CAF,  at 
delay  t  «  (2q  +  1)a;  again,  referring  to  (34), 

xa(y,(2q  +  1)a)  =  A^  exp[-12iruA(m  *  ^)]  Rfl  ((m  +•  |)A,(2q  f  1)a)  - 

m 

•  Jdt  exp(-12irvt)  R  (t,  (2q  +  1)a)  A  l(t  -(m  *  j)a)  « 

m 

-  x(v,(2q  *  1)a)  £  4<w  “  ?)  • 

m 

-£(-l)m  x(v  -  ®,  (2q  +  1)a)  for  all  w  .  (57) 

m 
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The  t  increment  is  again  a,,  but  it  is  shifted  by  a/2,  in  keeping  with  (39b) 
and  figure  4.  We  are  now  taking  a  horizontal  slice  at  t  -  (2q  f  l)a  in 
figure  4.  The  (-l)m  factor  Is  explained  by  (44). 

The  aliasing  lobes  In  (56)  and  (57)  are  separated  only  by  au  *  1/a  and 
will  overlap  on  the  v  axis  unless  a  <  (2F)”1;  see  (38)  and  figure  7.  Thus, 
the  approximate  CAP,  x,(v.na)  for  n  Integer,  suffers  overlap  due  to 

a 

aliasing,  just  as  the  WDF ,  W,(5-  a,  f)  for  n  Integer,  does;  elimination 
of  overlap  Is  achieved  only  If  the  stringent  requirement  a  <  (2F)^  Is  met. 
Furthermore,  again,  we  cannot  directly  average  the  two  results  In  figure  7, 
In  hopes  of  canceling  the  close-in  lobes  centered  at  «  •  t  1/a,  because  the 
two  delays,  2qa  and  (2q  1)6,  are  not  Identical. 

Equations  (56)  and  (57),  together,  Illustrate  that  approximations  to  the 

continuous  CAF  are  available  at  discrete  delay  values  with  separation  At  «  a 

and  at  a  continuum  of  v  values.  Now  we  have  succeeded  In  eliminating  the 

discrete  nature  of  the  other  of  the  two  Initial  time  variables  In  the  TCF, 

namely  t.  The  remaining  Fourier  transform  Into  the  SCF  domain  will 
eliminate  the  other  discrete  variable. 
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Figure  7.  CAF  Appro 
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APPROX  IMA  I  h  SPECTRAL  CORRELATION  f-UNCTION  VIA  CAP 


Guided  by  line  2  of  definition  (35),  we  obtain  the  following 
approximation  to  the  original  SCF  A(v,f): 

A^(w,f)  3  a  ^  exp( -I2irfna)  xa(v,nA)  * 


n 

exp(-12irfnA)  xa(v.nA) 

Q 


for  all  o,f  . 


(58) 


The  superscript  c  on  approximation  A,  denotes  that  we  are  obtaining  this 
result  by  way  of  the  CAF.  The  Increment  on  t  In  (58)  Is  a  -  A,  In 

T 

keeping  with  the  available  CAF  values  In  (56)  and  (57)  together. 


Let  n  -  2q  In  the  first  sum  In  (58),  and  let  n  »  2q  +■  1  In  the  second 
sum.  There  follows,  upon  use  of  (56)  and  (57), 

A^C)(w,f)  -  A  ^  exp(-i2irfA2q)  xa(v,2qA)  + 

q 


♦  a  2  exp[-12irf A( 2q  +1)]  xa(v,(2q  +  1)a)  - 

q 


^  exp(-12rrf A2q)  x(v  -  J,  2qa)  * 
m  q 


+  £(-1)"'  exp[  12irfA(2q  +  1)]  x  (v  -  J, 
m  q 


(2q  *  DA 


)■ 


(59) 
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But 

A  exp(  -I2irf  A2q)  x( » 1  »2qA)  «= 

q 

■  j* dt  exp( -1 2irf t )  x(v'»D  A  A(t  -  2qA)  ■ 

q 

-  A(w',f)®|  ^>4(f  -  gj)  - 

q 

-  \  ^A(v',f  -  j}')  ,  (60) 

q 

while 

a]>  exp[-l2irf  A(2q  +  1)]  x  (w* ,  (2q  +  da)  - 

q 

■  J dT  exp( -1 2-irf t )  x(u'.t)  a]!>a(t  -  (2q  f  Da)  ■ 

q 

-  A(v\f)  4^(“1)q  6(f  ‘  2A>  " 

q 

-  £2(“1)q  A(w',f  “  2A>  ’  (61)- 

q 

The  use  of  these  two  relations  In  (59)  yields  the  approximate  SCF 
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A<c)(»,n  -  \  ]>  £  A<”  ■  A-  f  -iJ>  - 

m  q 

♦i 

m  q 

-  21  2  A<«  -  J.  f  -  2/t>  for  all  *.f  •  (62) 

m  q 
m+q  even 

But  this  result  is  Identical  to  the  approximate  SCF  A^(»,f)  given 
in  (51)  and  figure  6.  That  Is,  we  obtain  the  same  continuous  approximation 
In  the  v,f  domain,  whether  we  approach  It  via  the  WDF  or  the  CAF.  This 
apparently  fortuitous  result  Is  due  to  the  fact  that  we  used  alj,  the 
available  Information  about  the  TCF  when  we  started  with  (39),  and  kept  all 
of  It  In  passage  through  the  WOF  or  CAF  domains. 

Figure  6  Is  again  applicable,  and  we  now  see  that  we  can  drop 
superscripts  w  and  c  from  (51)  and  (62),  respectively,  since  there  Is  only 
one  approximation  In  the  SCF  domain.  (The  comments  following  (51)  are  also 
directly  applicable  here.) 

A  rigorous  proof  of  the  equality  of  the  two  approximations  available  for 
the  SCF  Is  given  In  appendix  B.  It  utilizes  an  Impulsive  sampling  approach, 
similar  to  (9)  but  In  two  dimensions,  and  can  be  considered  as  an  alternative 
to  the  approximation  approach  developed  1 n  this  section.  Of  course,  the  end 
result  for  the  SCF  In  the  v,f  domain  is  again  (51)  or  (62)  or  figure  6. 


33 


TR  8533 


SUMMARY  STATUS  IN  ALL  FOUR  OOMAINS 

The  results  for  the  approximations  to  the  TCF,  WDF,  CAF,  and  SCF  are 
sketched  In  figure  8.  These  plots  are  a  condensation  of  the  exact 
analytical  results  given  by  (39),  (42)  l  (43),  (56)  £(57),  and  (51)  <£  (62), 
respectively.  For  example,  the  approximate  WDF  In  the  lower  left  of  figure  B 
Is  available  only  alon*  tfce  slices  where  t  *=  na/2,  n  Integer.  Along  these 
slices,  the  a  1  las  1  .*  lobes  (in  frequency)  alternate  In  polarity  if  n  Is  odd, 
but  remain  positive  for  n  even.  (Positive  lobes  are  drawn  toward  the  right 
side  In  the  figure). 

Horizontal  movement  from  one  diagram  to  another  In  figure  8  Is 
accomplished  by  a  Fourier  transform  from  t  to  v  (or  vice  versa). 

Vertical  movement  Is  according  to  a  Fourier  pair  relating  t  and  f. 

Finally,  the  diagonal  connection  between  R  and  A  ,  or  between  Wa  and 
x.i  Is  by  means  of  a  double  Fourier  transform. 

u 
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RECOVERY  OF  ORIGINAL  CONTINUOUS  TWO-DIMENSIONAL  FUNCTIONS 

We  have  seen,  by  means  of  (52) -(54)  and  figure  6,  that  the  original  SCF 
can  be  recovered  from  the  approximate  SCF 

Aa(u,f)  *  ^  ^  A(u  '  !•  f  “  for  a11  v’f  •  (M> 

q  m 
q+m  even 

by  means  of  the  diamond  gating  function  D(v,f)  in  (53),  provided  that 
a  <  T/F,  We  have  used  (51)  and  (62)  here,  and  dropped  the  superscripts  In 
accordance  with  the  discussion  following  (62).  This  means  that  we  have  the 
possibility  of  evaluating  the  original  continuous  TCF  R( t f t ) ,  WDF  W(t,f), 
and  CAP  x(v.t)  at  any  argument  values  we  please. 

SIMPLIFICATION  OF  SCF  Aa(v,f) 

It  would  be  an  extremely  tedious  task  to  evaluate  the  approximate  SCF 
A  directly  by  its  definition  (46)  coupled  with  (40)  and  (41),  which,  in 

A 

turn,  are  based  upon  starting  information  (3S).  In  fact,  there  Is  a 
startlngly  simple  way  of  computing  A_. 

Recall  from  (16)  that  spectrum 

S(f)  »  a  ^  exp(-12wfak)  s(ka)  - 
k 

S(f  -  £)  for  all  f  .  (64) 

k 
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Therefore 

5<f  *■  f)  S*(f  -  J)  -'2Xw  *  f  -  *)  s*(f  -*-}>. 

k  l 

A(V  ‘  *f  “  i0r  411  W,f  *  (65) 

k  i 

where  we  used  (35).  Now  let  m  *  k  -  l  and  q  -  k  f  l;  then 

rc ■  ~ -&  m  k,  ^""2  ^  ■  t  ,  meaning  that  m  ±  q  must  always  be  even.  Therefore, 
(65)  can  be  expressed  as 

S(f  ♦  J)  5*(f  -  f)  “  ^  ^  A<w  “  ?•  f  '  for  a11  w,f  ‘  <66> 

q  m 
q+m  even 

But  (66)  Is  Identical  with  (63)!  Thus  we  have  the  compact  result  for  the 
approximate  SCF 

Aft(v,f)  -  3(f  +■  {)  5*(f  -  f)  for  all  v,f  ,  (67) 

where 

5(f)  -  exp(-12irfak)  s(kA)  for  all  f  ,  (68) 

k 

in  terms  of  the  original  time  samples  {s(kA)|  . 


It  is  convenient  at  this  point  to  define,  for  all  f,  the  function 


S(f)  -  5(f)  rect(af)  - 


*  * 
exp(-12trfAk)  s(kA)  for  |f|  <  gj 


0  otherwise 


,  (69) 
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which  can  be  computed  directly  from  the  samples  {s(ka)}  .  Then  since 
a  <  l/F,  reference  to  figure  2  and  (18)  reveals  that 

5(f)  «  S(f)  for  all  f  .  (70) 

The  only  reason  for  distinguishing  between  T  and  S  Is  that  we  think  of  S  as 
being  computed  directly  from  samples  {s(kaj}  via  (69),  whereas  we  think  of  S 
as  being  computed  from  s(t)  via  Fourier  transform  (13).  Strictly,  since  S(f) 
Is  bandllmlted  to  t  F/2,  5"(f)  In  (69)  only  needs  to  be  computed  In  that 
somewhat  smaller  range  of  f. 

At  this  point,  we  refer  back  to  ( 52) -( 54)  and  figure  6  to  find  that 

A(v,f )  -  Aa(u,f)  0(v,f)  -  S(f  +  f)  S*(f  -  J)  for  all  v,f  ,  (71) 

since  only  the  origin  lobe  In  figure  6  can  contribute,  and  there  Is  no 
overlap.  Thus  we  have  a  very  direct  way  of  recovering  the  original  SCF  from 
the  time  data  samples:  compute  T(f)  from  (69),  and  then  A(«,f)  from  (71). 

All  these  results  are  predicated  on  sampling  rate  condition  a  <  l/F;  they  do 
not  require  A  <  (2F)"1 . 

If  we  substitute  (70)  In  (71),  we  hive  original  definition  (35).  Thus  we 
have  come  full  circle  on  the  SCF,  retuning  with  an  obvious  relation.  This 
Indicates  that  a  shortcut  could  have  bn»n  taken  with  regard  to  getting  the 
key  result  (67).  We  have  pursued  the  longer  route  because  It  Indicates  what 
the  complete  set  of  fundamental  processing  equations  are,  and  It  clarifies  a 
number  of  points  that  have  been  under  contention  In  the  literature. 
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RECOVERY  Of-  ORIGINAL  WDE  AND  CAE 

From  (37)  and  (71),  we  obtain  the  original  continuous  WDF  as 

W(t.f)  -  Jdv  exp(12irut)  3(f  +  f )  §*(f  -  |)  for  all  t.f  ,  (72) 

where  3(f)  Is  given  by  (69)  In  terms  of  samples  {s(ka)}  .  The  truncation  of 
3(f)  at  f  -  ±  (2a)"1  in  (69)  is  what  prevents  all  the  distant  si  del obes  of 
A4 ( v , f )  from  contributing.  We  could  hardly  have  expected  a  simpler  result. 

From  (37)  and  (71),  we  also  obtain  the  original  CAF  according  to 
x(v.t)  «  Jdf  exp(  1 2irf t )  S(f  ♦  |)  S*(f  -  f)  for  all  v.t  .  (73) 

Thus,  both  the  WDF  and  the  CAF  can  be  recovered  by  single  Fourier  transforms 
of  the  same  product  function,  but  on  complementary  variables  v  and  f, 
respectively. 

RECOVERY  OF  ORIGINAL  TCF 


Probably  the  best  way  to  recover  the  original  TCF  R(t,x)  is  by  means 
of  a  combination  of  (32),  (5),  and  (70): 


with 


R(t,r)  -  S(t  +  f)  s  (t  -  J)  , 


s(t)  -  Jdf  exp(12irft)  S(f)  . 


(74) 

(75) 


All  the  above  procedures  employ  S(f)  and  a  Fourier  transform  In  some 
fashion.  The  quantity  3(f)  can  be  computed  at  any  f  of  Interest,  directly 
from  samples  {s(kftj} ,  by  means  of  (69). 
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DLRE.C1  TIME  DOMAIN  RECOVERY  OF  CONTINUOUS  WDF 

We  have  given  two  alternatives  for  the  recovery  of  continuous  time 
waveform  s(t)  from  samples  {s(kfl)}.  They  are  (21)  or  (75)^(69).  If  we 
employ  the  former  In  the  time  definition  of  the  original  WDF  (line  1  of 
(26) ) ,  we  find,  for  all  t,f , 

W(t,f)  ■  J  dt  exp(  -1 2irf t)  s(t  +  |)  s*(t  -  f)  “ 

■»  a  £  £  expH2irfa(k  -  i)]  s(ka)  s*(u)  w0(t  -  m)  ,  (76) 

k  i 

where 


This  result  Is  equivalent  to  [5,  (5)4(&)3*  However,  as  noted  there,  this 
alternative  for  the  WDF  Is  not  computationally  attractive,  although  (76)  Is 
certainly  alias-free  because  It  Is  restoring  W(t,f)  Itself,  and  not  some 
approximation  to  It. 

As  an  aside,  If  the  frequency  domain  version  of  the  WDF  Is  used  Instead 
(line  2  of  (26)),  and  If  (21)  Is  Immediately  transformed  Into  the  frequency 
domain,  we  get  directly 
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S(f)  *=  3(f)  rect(Af)  =  S(f)  ,  (78) 

In  complete  agreement  with  (72). 

In  the  sequel  to  (62),  it  was  mentioned  that  an  alternative  approach 
Involving  Impulsive  sampling  could  be  used  to  get  various  Impulsive 
two-dimensional  functions  from  samples  {s(kaj} .  In  a  similar  vein,  the 
continuous  two-dimensional  functions  can  be  recovered  by  direct  convolution 
(Interpolation)  in  the  domaln(s)  of  Interest.  These  alternative  forms  are 
not  as  numerically  useful  as  the  ones  presented  above,  and  so  are  deferred 
to  Appendix  C.  However,  some  useful  insight  Into  the  Inadequacy  of  some 
past  attempts  at  Interpolation  Is  gained  by  this  alternative  viewpoint,  and 
the  readers  attention  is  directed  to  those  results. 

DISCUSSION 

In  retrospect,  (72)  and  (69)  are  an  obvious  result.  We  know  that  the 
original  continuous  WDF  Is  given'  by  (line  2  of  (26)) 

W(t,f)  -  j*  dw  exp(  12irvt)  S(f  +  |)  S*(f  -  jf)  for  all  t,f  .  (79) 

So  If  we  can  get  S(f)  exactly  from  samples  [s(kA)}  ,  In  some  (anv)  fashion, 

we  can  get  W  via  (79).  But,  In  fact,  3(f)  In  (69)  Is  Identically  equal  to 

-1 

S(f )  for  all  f,  when  a  <  1/F.  Condition  a  <  (2F)  Is  patently  unnecessary 
and  too  restrictive.  A  similar  comment  holds  with  regard  to  the  CAF. 
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Given  samples  £s(kA)},  the  function  5(f)  In  (69)  can  be  computed  at  any 
leslred  f  values  of  Interest.  Therefore  the  product 

S(f  +  })  f(f  -  })  (80) 

required  In  (72)  and  (73)  can  be  computed  at  any  v,f  values  needed,  and  the 
Integrals  for  W(t,f)  and  x(v*t)  evaluated  very  accurately  at  any 
arguments  of  Interest. 

This  Is  the  major  difference  relative  to  the  TCF  R(t,x),  (74),  which 
could  only  be  calculated  at  Interspersed  points  In  the  t,T  plane  from  the 
available  data;  see  (39)  and  figure  4.  Strictly,  waveform  s(t)  could  be 
Interpolated,  and  then  TCF  R(t,t)  filled  In  at  the  Intermediate  points  of 
Interest  In  figure  4.  This  viable  alternative  requires  just  slightly  more 
calculations  that  the  frequency  domain  approach  given  above;  we  will  discuss 
and  compare  both  alternatives  In  a  later  section. 

In  practice,  5(f)  will  only  be  calculated  at  a  discrete  set  of 
frequencies,  in  order  to  economize  on  computational  effort.  We  then  find 
that  the  product  function  (80)  Is  available  at  Interspersed  points  In  the 
v,f  plane  In  an  Identical  manner  to  that  for  TCF  R(t,r)  In  figure  4. 

In  fact,  If  5(f)  Is  computed  only  for  f/af  ■  Integer,  then  (80)  Is 
available  only  at 

f/af  -  Integer  ,  v/af  ■  even  Integer  (81) 

and  at 

f/af  *  odd  lnteger/2  ,  v/af  ■  odd  integer  .  (82) 

Just  as  this  type  of  Interspersed  sampling  required  a  finer  sampling  Interval 
In  the  time  domain  (see  figure  5  and  (45)),  so  also  Is  a  finer  Increment 
required  here  In  the  frequency  domain.  Namely,  we  must  have 
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Af  <  21  •  (83) 

in  order  to  avoid  aliasing  In  the  t  domain  of  the  reconstructed  WDF  via 
Fourier  transform  (72).  The  same  requirement  holds  for  aliasing  control  In 
the  t  domain  of  the  reconstructed  CAF  via  (73).  Here,  T  Is  the  overall 
effective  duration  of  waveform  s(t): 

|s(t)|  a  0  for  jt|  >  T/2  .  (84) 

See  figure  9.  (The  waveform  can  be  centered  at  t  «  0  without  loss  of 
generality,  merely  by  time  delaying  It.)  However,  there  Is  a  very  convenient 
and  efficient  way  to  meet  requirement  (83),  as  will  be  shown  shortly,  whereas 
requirement  (45),  a  <  (2F)"1,  is  very  unattractive,  at  least  through  direct 
sampling  of  time  domain  waveform  s(t). 

Since  the  total  extent  of  the  spectrum  S(f)  Is  F  Hz  (see  (14)  and 
figure  1),  waveform  s(t)  cannot  be  strictly  time-limited.  However,  we 
assume  that  a  finite  T  value  can  be  found  for  figure  9  such  that 
approximation  (84)  is  a  good  one.  Strictly,  (84)  should  read 

- «  1  for  Jt|  >  T/2  .  (85) 

max  }s(t)] 
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DISCRETE  PROCESSING  IN  FREQUENCY  DOMAIN 


Thus  far,  we  have  not  discretized  In  the  frequency  domain;  both  v  and  f 
have  been  allowed  to  take  on  continuous  values.  So  we  have,  from  (69), 


§(f)  -  a]>  exp(-12irfak)  s(ka)  for  |f|  <  ^  .  (86) 

k 

Now,  suppose  that  we  only  evaluate  S(f)  at  a  set  of  discrete  frequencies, 
according  to 

S(jjJ)  -  a  ^  exp(-12irnk/N)  s(ka)  for  |n|  <  |  .  (87) 

k 


Since  the  sum  on  the  right-hand  side  of  (87)  has  period  N  In  n,  we  can 
evaluate  It  quickly  via  a  collapsed  N-poInt  FFT ;  see  ( 1 0) -( 1 2) .  The 
negative  n  values  desired  In  (87)  are  easily  accommodated  by  means  of  a 
modulo  N  look-up  In  the  FFT  output. 

EVALUATION  OF  WDF 


The  Increment  In  argument  f  of  3(f)  In  (87)  Is 

*f  ■  si  • 


(88) 


In  order  to  use  these  results  In  approximating  Integral  (72)  for  the  WDF, 


W(t,f)  -  Jdv  «xpM2««t)  5(f  +  J)  §*(f  -  £)  , 


(89) 


we  need  to  have  the  Increment  In  v  satisfy  (due  to  the  v/2  arguments) 
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1  1  2 
r\  A  S  ll.  |  1*6*|  A  II.  ' 

2  v  NA  u  NA 

So  If  we  limit  v  In  (89)  to  the  values  ,  one  possible  Trapezoidal 
approximation  to  W  Is 

W(t.f)  i  jjf  2  «*P<12*  NA  §(f  f  NA*  §*<f  -  HA>  for  a11  t*f  • 


(90) 


(91) 


& 

This  function  Is  defined  on  a  continuum  In  time, frequency  space.  But  this 
can  be  developed  according  to 


SJ(t.f)  -  [dv  exp(12irvt)  S(f  *  J)  §*(f  -  f)  jjf  2  6(v  “  nJ)  " 

& 

-  W(t,f)®5  4(t  -  l  %)  ■ 

l  4 

■  J>W(t  -  l  f)  for  all  t,f  ,  (92) 

i 


where  we  used  (72). 


Since  waveform  s(t)  Is  approximately  limited  to  (t|  <  T/2  (see  (84)  and 
figure  9),  then  WOF  W(t,f)  Is  also  approximately  limited  to  |t|  <  T/2,  as 
may  be  seen  from  line  1  of  (26).  The  approximation  (92)  then  appears  as  In 
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In  order  that  approximate  WOh  W  not  suffer  significant  overlap  In  time, 
we  see  that  we  need 


1 

2 


< 


M  _  I 
2  2  * 


i  .e. , 


(93) 


The  FFT  size,  N,  in  (87)  must  be  at  least  twice  as  large  as  the  number  of 
samples,  T/a,  taken  of  waveform  s(t).  Recalling  (88),  this  Inequality 
becomes 


a 


f 


J.  <  -1 

Na  2T 


(94) 


consistent  with  (83),  as  predicted.  Thus,  approximation  W(t,f)  Is  an 
extremely  good  approximation  to  W(t,f)  for  jt|  <  T/2  If  FFT  size  N  satisfies 
(93).  The  goodness  of  ft  depends  critically  on  the  degree  of  satisfaction  of 
(85). 


More  generally,  we  could  limit  the  »  values  In  (89)  to 

•  ■  »0  +  U  '  "o  ,rbUr4rV  <\  -  si>  .  <“) 

getting  alternative  approximation 

EJa(t,f)  =  Jdv  exp(l2irtv)  S(f  +  |)  §^( f  -  J)  (jf  ^  4(w  -  w0  “  Na^  " 

s. 

-  W(t,f)®  jexp^irtv^ajt  “  1  “ 

% 

W(t  “  1  **2*  f)  e«p(iifv0Nat)  for  all  t.f  .  (96) 

i 

The  plot  for  | W a ( t , f ) j  Is  Identical  to  that  of  jft(t,f)]  in  figure  10,  since 

the  magnitude  of  phase  factor  exp(1irt>oNat)  Is  1  for  all  l,  The  main  lobe, 

t  -  0,  Is  unaffected  by  the  choice  of  #  .  So  criterion  (93)  Is  again 

o 

sufficient  to  avoid  time-aliasing  In  tJ  ,  regardless  of  shift  v  . 

A  0 
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DISCRF.il/A110N  IN  1  IMF.  AND  FREQUENCY 


For  convenience,  we  therefore  return  to  W  in  (91)  and  get,  In  particular, 
the  values 


-  jji  2  #xp(12ir  &  t)  S(^)  s"(2~i)  for  .11  t  .  (97) 

l 

where  we  must  choose  frequency  f  -  ,  In  order  to  use  the  available 

samples  of  S  In  (87).  We  now  further  choose  time 


t  .  !M  N 

1  2  M 


and  get  the  approximation  in  the  form 


(98) 


».RI)  ■  jj£  2  0«PO2«»l'W>  .  (99) 

t 


The  reason  for  this  choice  of  t  values  is  that  this  sum  can  be  accomplished 
as  a  collapsed  M-poInt  FFT  as  described  in  (10)-(12).  The  range  of  values 
that  must  be  covered  is 


MA  U  <  I 

2  M  2 


i.e.,  |ml  <  |j| 


and 


JnL  <  E 

Na  2  ’ 


i.e.,  |n[ 


(100) 

(101) 


Coincidentally,  this  identical  procedure  above  has  already  been  derived 
by  the  author  in  [9,  (A-13)I(A-14) ].  However,  it  was  done,  at  that  time,  to 
generate  slices  in  time  of  the  WDF ,  without  realizing  that  the  procedure  a : so 
had  an  alias-elimination  feature.  Requirement  (93)  was  [9,  (A-5)£ (A-6)  j . 
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INCME.Mb.NTS  IN  t  AND  f 

The  increments  In  t  and  f  In  approximation  W(t,f)  in  (99)  are 

4t  ■  f  &  aml  4f  ■  si  ■ 

Since  the  original  WDF  is  given  by 

W(t,f)  -  Jdr  exp( ~i2irf t)  R(t,t)  , 

and  the  effective  extent  of  R(t,t)  In  t  Is  ±T  for  the  waveform  s(t) 
satisfying  (84),  then  we  must  require 

4f  <  2t  •  *•••»  N  >  21  ,  (104) 

in  order  to  track  the  variation  of  W(t,f)  in  f,  (See  appendix  A.)  However, 
this  was  a  condition  already  encountered  in  (93). 

Furthermore,  since  we  have  the  alternative  Fourier  transform 

W(t,f)  -  J  dv  exp(12irvt)  A ( u , f )  ,  (105) 

and  the  extent  of  A(v,f)  in  v  Is  ±F.  then  we  must  also  have 

4^  ^  2^  «  1.e.,4<jrjJ,  (106) 

in  order  to  track  W(t,f)  variations  in  t.  Now  If  we  choose  M  smaller  than  N. 
say  M  ■  N/2,  then  we  obtain  condition  a  <  (2F)"\  Put  this  Is  a  finer  time 
sampling  Increment  than  required.  Also  (102)  gives  at  »  a,  which  does  not 
track  W(t,f)  adequately  in  time;  see  ( 24 ) -( 31 ) .  Conversely,  If  we  choose  M 
larger  than  N,  say  M  •  2N,  then  we  get  a  <  2/F,  which  is  already  accommodated 

by  the  earlier  requirement 
48 


(102) 


(103) 
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A  <  1/F.  And  -•  a/4,  which  is  overly  fine  in  time.  So,  In  order  to 
minimize  the  range  of  m  values  needed  for  investigation,  we  choose 

M  -  N  .  (107) 

Notice  that  the  time  Increment  At  In  (102)  for  fl(t,f )  Is  then  At  «  A/2, 
not  the  a  that  was  sufficient  for  sampling  s(t).  This  is  consistent  with 
the  fact  that  W(t,f)  can  be  sharper  in  t  than  s(t). 

SUMMARY  OF  WQF  EQUATIONS 

Here  we  summarize  the  major  assumptions  and  requirements  and  list  the 
major  equations  by  which  the  approximate  WDF  W  can  be  computed. 


Assumptions: 

S(f)  -  0  for  If)  >  F/2  , 

|s(t)|  %  0  for  |t|  >  T/2  .  (108) 

This  Is  no  loss  of  generality,  since  the  waveform  can  be  time  delayed  and 
frequency  shifted  as  desired. 

Requirements: 

(109) 

(110) 

49 


Equations: 


5i> 


N  >  &  (>  2TF)  . 

A 


A 2.  exp(-12imk/N)  s(kA) 
k 


0  otherwise 


for  )n|  <  £ 
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W( 


034  _4. 

2W 


iii  Z  e^PC^^iL/N)  Sl^)  S*!0-^1) 

i 

for  lm\  <  J  ,  |n|  <  |  Fa  . 


(Ill) 


Operation  (110)  can  be  accomplished  by  a  single  N-poInt  FFT  with 
collapsing,  while  (111)  requires  an  N-poInt  FFT  for  each  n  value  of  Interest. 
The  latter  N-poInt  FFT  (for  a  given  n)  will  sweep  out  N  values  of  integer  m} 

this  will  cover  a  total  time  range  of  N  £  >  T  ,  as  desired.  The  total  number 
of  n  (frequency)  values  to  be  searched  Is  NFd.  Values  of  ft  for  negative 
values  of  m  are  available  In  location  m  modulo  N  of  the  FFT  output. 


For  the  most  advantageous  choices  of 

A  -  1  ,  N  «  «  2TF  ,  (112) 

we  have  ranges 

|mj  <  TF  ,  |n l  <  TF  . 


EVALUATION  OF  CAF 


The  equations  for  calculation  of  the  CAF,  from  discretized  frequency 
samples  of  5(f),  are  very  similar  to  those  given  above  for  the  WDF. 
Accordingly,  they  are  deferred  to  appendix  0. 
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INTERPOLATION  OF  TIME  WAVEFORM 

It  was  mentioned  In  the  sequel  to  (80)  that  available  data  samples 
{s(kA)1  could  be  Interpolated  In  time,  to  fill  In  the  vacant  spots  of  TCF 
R(t,r)  In  figure  4.  One  procedure  to  accomplish  this  Is  by  direct  use  of 
sinc-function  Interpolation  in  (21).  However,  a  more  efficient  procedure  Is 
by  use  of  FFTs. 

From  (22),  we  can  get'  samples  of  the  approximate  spectrum  according  to 

S(jjJJ)  -  exp(-12irmk/M)  $(ka)  for  |m|  <  |  ,  (113) 

k 

which  can  be  accomplished  by  an  M-poInt  FFT.  The  frequency  increment  Is 
af  ■  (Ma)"1.  The  original  continuous  time  waveform  Is 

s(t)  -  Jdf  exp(12irft)  S(f)  - 

(2A)""1 

df  axp(12irft)  5(f)  » 

M/2 

:M4  ^  •xp(l2rS“t)  5(jj*)  ■ 

m— M/2 

5  s(t)  for  all  t  ,  (114) 
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where  we  have  utilized  the  fact  that  a  <  l/'F.  (The  tic  mark  on  the  summation 
symbol  Indicates  that  the  summand  values  at  m  ■  tM/2  must  be  scaled  by  l/2#) 
ETut  approximation  $(t)  can  be  developed  In  the  form 

J. 

20  M/2 

S(t)  -  J  df  exp(12nft)  5(f)  ^  ^  6(f  -  Jjj)  - 

zL  m-572 

2a 

-  j'df  txp(12«ft)  S(f)  <(f  '  Sa>  ■ 

m 

s(t)®^  A(t  -  mMa)  s(t  -  mMa)  for  all  t  .  (115) 

m  m 

This  aliased  function  Is  Illustrated  In  figure  11. 


i(B 
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In  order  to  avoid  time  overlap  in  figure  11,  we  must  have 

M  >  J  (>  TF)  .  (116) 

A 

If  this  FF r  size  requirement  for  (113)  Is  met,  then  figure  11  reveals  that 

$( t)  -  s(t)  for  |t|  <  ^  ;  (117) 

that  Is,  we  can  expect  that  l(t)  Is  a  good  approximation  to  s(t). 


In  particular,  If  we  want  to  Interpolate  samples  {$(kA)}  by  a  factor  of 
2,  we  have,  from  (117)  and  (114), 


M/2 


s(^)  =  S(^)  -  ^  exp(12ir  g)  §(jg)  for  Jn|  <  M  .  (118) 


m— M/2 


This  can  be  done  by  a  2M-po1nt  FFT  with  rero-flll,  and  made 

available  by  the  M-poInt  FFT  In  (113),  The  two  requirements  that 
must  be  met  are 


<  1 
F f 


4  -  c* 


2M  >  (>  2TF)  . 


(119) 


This  Is  the  procedure  utilized  In  [8], 
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DISCUSSION/SUMMARY 


We  have  two  alternatives  for  computing  an  alias-free  WDF  and  CAP,  The 
requirements  that  must  be  met  are 

A  <  p  ,  N  >  |  ,  N  •  2M  >  ^  .  (120) 

The  processing  equations  are  summarized  below, 


TIME  DOMAIN  APPROACH 


exp(-12irmk/M)  s(kA)  for  |m|  <  |  , 


(121) 


M/2 


Rp  •*P(12trmn/N)  S(rJ)  for  |n|  <  M 

m— M/2 


M) 


0  oth«rwl„ 


(12?) 


^(®2 «  fjJ)  -  A  ^  exp(-l2irni/N)  A)  A) 


for  |m|  <  J  ,  |n|  <  |  FA  <  ^ 


2  1 


(123) 


In  the  last  equation,  we  obtain  slices  of  the  WOP  In  frequency  (n)  at  fixed 
time  (m).  (This  result  Is  In  essential  agreement  with  [7,  (77)  and  (74)]; 
however,  we  are  not  restricted  here  to  T  ■  MA.) 
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FREQUENCY  DOMAIN  APPROACH 


(124) 


^(!D2*  Na5  “  N 1%  «xp(  12wml/N)  5  (IL^1) 

l 

for  |mj  <  J  <  |  ,  )n|  <  J  FA  .  (125) 

Here,  we  obtain  slices  of  the  WOF  In  time  (m)  at  fixed  frequency  (n),  (This 
result  Is  In  essential  agreement  with  [7,  (81)];  however,  we  are  not 
restricted  here  to  2T  -  NA.) 


COMPARISON 

Since  waveform  s(t)  -  0  for  [tl  >  T/2,  the  sum  on  k  In  (121)  can  be 
confined  to  )k)  <  2a  <  2  "  4  *  the  sum  on  k  4|n  O24)  can  &• 

confined  to  |k|  <  ^  c  |  . 

The  time  domain  approach  requires  cne  M-polrit  (M  *  N/2)  F FI ,  one  N-poInt 
FFT ,  and  then  an  H-poInt  FFT  for  each  time  Index  m  of  Interest.  Thr 
frequency  domain  approach  requires  one  N-poInt  FFT  and  then  an  N-poInt  FFl 
for  each  frequency  Index  n  of  Interest.  Thus,  the  time  domain  approach 

requires  one  additional  M  ■  ^  -  point  FFT,  which  Is  a  negligible  difference, 
compared  with  the  multiple  FFTs  that  must  be  performed  for  (123)  or  (125). 
The  major  difference  between  the  two  approaches  appears  to  be  whether  one 
wants  slices  In  frequency,  or  slices  In  time,  of  the  WDF. 
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It  should  be  emphasized  that  the  procedure  here  reconstructs  the 
original  WDF  of  the  continuous  waveform  s(t),  at  any  time, frequency  (t,f) 
arguments  of  Interest,  from  the  samples  [s(ka)j  .  There  Is  no  need  to  define 
or  set  up  some  arbitrary  discrete  form  of  the  WDF.  The  discretization  of 
the  t,f  arguments  of  the  WQFs  Is  undertaken  only  after  this  reconstruction 
procedure  has  been  delineated  (via  the  v,f  plane)  and  the  sufficiency  of 
sampling  requirement  &  <  1/F  established.  Of  course,  this  evontual 
discretization  In  time/frequency  Is  necessary  In  order  to  reduce  the  general 
procedure  to  a  practical  efficient  algorithm.  A  similar  set  of  arguments 
applies  equally  well  to  the  reconstruction  of  the  CAF . 

An  alternative  philosophy,  to  developing  Trapezoidal  approximations  for 
the  various  two-dimensional  functions,  Is  given  In  appendix  B  In  terms  of  a 
pair  of  Interspersed  Impulsive  trains  In  the  t,T  plane.  The  end  result 
for  the  approximate  SCF  In  the  «,f  plane  Is  shown  to  be  Identical  to  that 
obtained  earlier. 

The  connection  of  this  two-dimensional  Impulse  train  with  direct  time 
domain  sampling  of  the  waveform  s(t)  Is  considered  In  appendix  E.  Again, 
the  two  approaches,  time  domain  sampling  versus  t,T  domain  sampling,  are 
shown  to  yield  the  same  result.  Finally,  the  fundamental  rules  and  patterns 
relating  two-dimensional  interspersed  Infinite  sampling  trains  are  displayed 
at  the  end  of  appendix  E. 
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APPENDIX  A.  EXTENTS  AND  PATES  OP  VARIATION 
OF  TCP.  WDF ,  CAF,  SCF 

From  (5),  we  have  Fourier  pair 

S(f)  -  j' dt  exp(-12irft)  s(t)  , 

s ( t )  -  y df  exp(l2irft)  S(f)  {  (A-l ) 

and  from  figures  1  and  2,  we  know  that  If 

S(f)  -  0  for  If i  >  F/2  ,  (A-2) 

then  samples  (s(kAt)}  are  sufficient  for  reconstruction  of  s(t)  If  at  <  1/F. 
That  is,  If  a  spectrum  Is  bandl 1ml ted  In  the  frequency  domain  to  total 
extent  F  Hz,  samples  of  the  corresponding  time  function  must  be  taken  with 
time  Increment  &t  <  1/F,  In  order  not  to  lose  any  significant  Information. 

In  a  similar  vein,  the  duality  of  the  equations  In  (A-l)  Indicates  that 
If,  Instead, 

s(t)  -  0  for  Jtl  >  T/2  ,  (A-3) 

then  spectrum  samples  {s(nbf)]are  sufficient  for  reconstruction  of  S(f)  If 
frequency  Increment  <  1/T. 

The  general  rule,  here,  Is  that  If  a  function  In  one  domain  Is 
essentially  limited  to  overall  extent  E,  samples  In  Its  Fourier  transform 
domain  must  be  taken  finer  then  1/E,  In  order  not  to  lo*e  any  Information. 
This  rule  will  be  used  frequently  below. 


A-l 


TR  8533 


EXTENTS  OF  TCF ,  WDF,  CAF,  AND  SCF 


Henceforth,  we  assume  that 

$(f)«  0  for  |f|  >  F/2 
and 

i(t)*  0  for  |t|  >  T/2  .  (A-4) 


Thus,  the  overall  frequency  and  time  extents  art  approximately  F  Hz  and  T 
seconds,  respectively.  It  then  readily  follows  from  (32)-(38),  namely, 


R(t,T) 

■  s(t  +  f)  »  (t  -  J)  , 

(A-5) 

A(vif ) 

-  S(f  +f)  S*(f  -  J)  , 

(A-6) 

W(t.,f)  • 

dT  exp(-12irfO  R(t,T)  * 

(A-7) 

• 

< 

J* dv  exp(12irvt)  A(v,f )  , 

(A-8) 

x(w.t)  ■ 

j*dt  exp(-12irvt)  R(t ,t)  » 

(A-9) 

■ 

j'df  exp< 1 2wf r )  A(v,f)  , 

( A-l  0 ) 

that  the  extents  of  these  functions  are  as  depicted  In  figure  A-l. 

The  solid 

curves  depict  the  contour  level  within  which  the  function  Is  essentially 

t2 

concentrated.  In  fact,  for  Gaussian  waveform  s(t)  -  a  exp(-  -k«),  the 

2  a 

choices  T  «  4a,  F  •  ,  for  example,  give  these  exact  results  In  figure  A-l, 

ita 


at  the  exp(-4)  -  .018  level. 


A-2 


Figure  A-l .  Extents 
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RAT tS  OP  VARIATION 


We  now  combine  the  sampling  rule  deduced  under  (A-3)  with  the  extents  In 
figure  A-l.  From  (A-7)  and  (A-8),  we  must  have  sampling  Increments  (l.e., 
spaclngs  of  argument  values  on  left-hand  side)  satisfy 


4f  <  Jt 


*  7F 


to  adequately  track  WOF  . 


(A-l  1 ) 


From  ( A-9)  and  (A-10),  we  need 

L 


t  <  i 

■t  F 


to  track  CAF  . 


By  Inverting  (A-7)-(A-10) ,  we  express 


R(t.-r)  -  j*(1f  exp(l2irft)  W(t,f)  - 
■  j*  dv  exp(12*vt)  w ,  t )  , 

A(v,f )  ■  ^dt  exp(-12irvt)  W(t,f)  • 


JdT 


exp(  -12irft)  x ( v •  t )  . 


(A~12) 


(A-13) 

(A  -14) 

(A -15) 

(A-l  6) 


Therw  follows  from  (A-l 3)  and  (A-l 4),  that  we  need 


\<f 


A  <  -1 

<  2F 


A -4 


to  track  TCF  , 


(A-l  7) 
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while  from  (A-15)  and  ( A  - 1 6 ) , 


A 


v 


A 


f 


to  track  SCF  . 


(A -18) 


These  last  four  restrictions  are  Identical  to  those  given  In  (A-11)  and 

(A-12).  The  total  number  of  samples  required  to  completely  describe  any  one 

2  2 

of  the  four  two-dimensional  characterizations  Is  TP. 


SAMPLING  RATE  FOR  REAL  WAVEFORM 

Requirement  (17)  on  the  time  sampling  Increment,  A  <  1/F,  for  recovery  of 
the  time  waveform  $(t),  Is  based  upon  figures  1  and  2  for  a  complex  envelope 
waveform  s(t).  if  s(t)  were,  Instead,  a  real  waveform  ( t) ,  the  earlier 
development  covers  this  case  as  well,  but  with  a  change  In  notation.  The 
spectrum  S^f)  of  waveform  s]  ( t)  Is  symmetric  about  f  -  0,  as  depicted  In 
figure  A-2.  F  Is  now  the  total  frequency  extent  of  the  positive-frequency 
components  of  s^ (t) . 

We  now  have  frequency  limit 

^1  P 

2  "  fc  +  2  •  ( A— 1 9 ) 

and  the  stringent  requirement  (45),  for  an  unaliased  WOF ,  becomes 

4  <  ^  -  *7777?  ■ 

For  a  narrowband  waveform,  f  »  F,  this  requires  an  unnecessarily  high 

V 

sampling  rate,  compared  to  what  would  be  required  for  the  waveform 


A-5 


corresponding  to  single-sided  bandwidth  F.  Extraction  of  the  complex 
envelope  (or  analytic  function)  of  s^t)  would  return  us  to  WDF  requirement 

a  <  jp  *  as  In  (45).  This  pre-processing  feature  Is  recommended  for  all  real 
waveforms.  However,  we  also  want  to  avoid  this  more  stringent  WDF 
requirement  and  be  subject  only  to  the  a  <  1/F  limitation. 


Figure  A-2.  Spectrum  of  Real  Waveform 


TR  8533 


APPENDIX  B.  IMPULSIVE  SAMPLING  APPROACH 


Instead  of  trying  to  develop  Trapezoidal  approximations  to  the  WDF,  CAF, 
and  SCF  Integrals  from  the  available  Information  about  the  TCF  In  (39),  we 
adopt  here  the  philosophy  that  continuous  TCF  R(t,r)  has  had  a  pair  of 
Impulsive  trains  applied  to  It,  yielding  the  Impulsive  approximation 


R^t.t)  5  A2  R(t,x) 


^  6 ( t  -  mA)^  *(t  -  2qA)  f 
.  m  q 


(m  4 • 


A*  R(t.r) 


n  L 
nH  even 


A(t  -  I1J)  4(t  -  IA)  - 


r 


R(flf,lA)  4(t 
n  l  £ 


n+L  even 


-  W)  «(T  -  LA)  . 


(B-l ) 


That  Is,  a  couple  of  two-dimensional  Impulse  trains,  Interspersed  In  the 
t.t  plane,  have  been  applied,  so  as  to  use  all  the  available  Information 
In  (39).  This  Is  identical  to  the  result  for  the  TCF  of  impulsive  time 
waveform  s^(t),  obtained  by  multiplying  s(t)  by  a  sampling  train;  see 
appendix  E. 


We  now  define  the  corresponding  WDF,  CAF,  and  SCF  as  rigorous  Fourier 
transforms  of  (8-1),  using  the  standard  forms,  for  all  argument  values, 


B-l 
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W^t.f)  -  y  dr  exp(  -12irf t)  R^t,*)  , 

Xi  (w,t)  -  Jdt  exp(-12irvt)  R.j  (t,t)  , 

A^(v,f)  -  Jdt  exp(-12*vt)  W^t.f)  ■ 

■  J*  dt  «xp( -I2nf t )  x^(v,t)  - 

dt  dt  exp(-12irvt  -  12irfx)  R^(t,t)  .  (B-2) 

These  exact  Interrelationships  Indicate  that  the  iftM  SCP  A,|  will  result 
from  TCP  R>,  whether  we  proceed  by  way  of  the  WDF  or  the  CAP. 

Me  have,  In  detail,  WDF 

W1  (t,f )  -  Jdt  exp(-12irft)  R^(t ,t)  - 

»  ^  a  A(t  -  mb) 
m 

+  |  a  a  (j  -  (m  +  |)a)  J  dt  exp(-12irft)  R(t,T)  2A^>  a(t  -  (2q  +  1)a)  ■ 
m  d 

f  ^ 

-  |  a]>  4(t  -  (DA)  W(t,f)®£i(f  -  jjj)  + 

m  Q 

+  £  A^  *(*  “  (»  +  |)a)  W(t,f)  ("1  )**  *(?  "  gj)  * 
m  d 

-  \  A  ^  *(t  -  IRA)  ^  W(t,f  -  j§)  + 
m  q 


I 


dt  exp(~12*ft)  R(t,t) 


26 1 
d 


S(t  -  2qA)  + 


+  2  » 


B-2 


S’  *(}  “  U"  +  J)A)£(-l>q  W(t,f  -  gA) 

m  q 


(B-3) 
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c  2  A(t  '  w<mA>f  “  2^)  * 

ID  q 

+  <(t  -  (m  +  J)a)£  H)q  W^m  +  *)4,f  -  gj)« 

m  q 

•  ^  4  ^j>  4(t  -  inA)  W#(m4,f)  +  ^  a  4  -  (m  +  Wa  +  ^)4,f)  * 

m  ffi 

-  *(t  -  “£>  Wa(Qf>f)  ,  (B-4) 

n 

using  (42)  and  (43).  Thus,  the  areas  of  the  Impulses  In  are  equal  to 
the  approximations  W,  developed  In  (42)  and  (43),  within  a  scale  factor  of 

A 

4/2. 

Continuing  on,  from  (B-3),  the  SCF  is 

A^(v,f)  ■  j* dt  exp(-12wvt)  W^(t,f)  * 

■  exp(-12«rvt)  W(t,f  “  2A>  4^  4(t  -  in4)  + 

q  m 

+  ££(-l)q  I  dt  exp(-12wvt)  W(t,f  “  2?)  4^  4^t  -  (m  f  l)a)  - 
q  m 


B-3 
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'  2  1  -  24>®2  4<”  ‘  !>  * 

q  m 

♦  |2<-i>q  A<»-f  ■  4(v  ■  !> ' 

q  m 

■  2  2  2  A(a  '  41  f  ‘  U>  * 
q  m 

*!I  £<-1>q*"  A(»  -B,  f  -  55)  - 

q  m 

-  y  2!  A(v  -  J,  f  -  2J)  for  all  v,f  ,  (8*5) 

q"  m 
q-Hn  even 

Thus,  the  SCF  Aj ,  resulting  from  the  Impulsive  sampling  approach  applied 
to  the  TCF ,  Is  not  Impulsive  at  all  In  the  v,f  plane,  and  Is  Identical  with 
the  approximations  A.  developed  In  (51)  and  (62). 

Q 


Proceeding  Instead  via  the  CAF,  we  have 


Xi(v.t)  ■  Jdt  exp(-12*vt)  R^(t,r) 


•  A  ^4(t  -  2qA)  J dt  exp(-12irvt)  R(t,t)  A  4(t  -  mA)  +• 
q  m 

A  2  *(*  "  (2<1  +  dt  txp(-12w*t)  R(t,t)  A^>  4  jt  -  (m  ■¥  j)a) 

a  m 


B-4 
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■  A  6(t  -  2qA)  x(v,t)®^  4(t>  -  ®)  ¥ 
q  in 

+  a£  *(r  -  <2q  ♦  1)A)  *(*,t)0£(-1)b  «(»  ~  J)  - 
q  m 

-  A^  4(t  -  2qA)^  x^w  “  ?•  f 
q  m 

♦  A^  a(t  -  (2q  +  Da)  x(*  *  }.  t)  -  (B-6) 

q  m 

“  A  4(  t  -  2qA)  x^»  “  2qA^  4* 

q  m 

♦  A  ^  |(t  -  (2q  +  1)^  £(-Dm  xfv  -  ®  (2q  +  1)a)  » 
q  m  '  ^ 

-  A  ^S(t  -  2qA)  xa(».2qA)  +  A^4  (r  -  (2q  4*  Da)  xa(y.(2q  4-  Da)  - 

q  q 


■  a)>  4(t  -  nA)  xa(».nA)  ,  (8-7) 

n 

via  (56)  and  (57).  Thus,  the  areas  of  the  Impulses  In  x^  are  equal  to 
the  approximations  x.  In  (56)  and  (57),  within  a  scale  factor  of  A. 

9 


B-5 
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Now  strictly  speaking,  there  Is  no  need  to  proceed  to  SCF  via  the 
CAF  Xj *  since  (B-2)  shows  that  there  Is  only  one  function,  regardless 
of  how  reached.  Nevertheless,  for  completeness,  we  also  present  the  last 
route.  We  have,  using  (B-6), 


Aj(v,f)  ■  J  dt  exp( "1 2 ir f t )  x^(u,t)  ■ 

2  ^  jdt  exp(-12irft)  x(«  -  ®.t)  2A^>6(t  -  2qa)  + 


+ 


dt  «xp( -1 2irf t)  x(v  -  J,t)  2A^T  4(t  -  (2q  f  1)a)  - 
m  ^  q 


*  2  £  A(w 
m 


♦il'-1*1" A<” 

m 


-  ID 

A' 


f)®2(-uq  «(f 
q 


■  f  -  si>  *!22<-i>'’+q  *«•-;.  f  -  j®>  - 

m  q  m  q 


-2  A(u  -  J,  f  -  5J)  for  all  u , f  .  (B-B) 

m  q 
m+q  even 


As  anticipated,  this  Is  Identical  with  (B-5).  Thus  we  get  a  unique  sen  In 
the  v.f  plane. 


B-6 
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APPENDIX  0.  RECOVERY  VIA  DIRECT  CONVOLUTION 

This  appendix  Is  closely  coupled  with  the  previous  one;  It  shows  how  to 
recover  the  original  continuous  two-dimensional  TCP,  CAF ,  and  WDF  from  their 
Impulsive  counterparts.  From  (B-5),  (B-8),  (51),  (52),  (53),  (54),  and 
figure  6,  the  original  SCF  Is 

A(w,f )  -  A^(v,f )  D(v.f)  .  (C-l) 

WDF  RECOVERY 

We  have,  using  (8-4) , 

W(t,f)  ■  J  dv  exp(12irut)  A^(w,f)  D(w,f)  ■ 

t 

-  W^t.f)  #d(t,f)  - 

*  ^  Wa(D2,f)  2  d(t  '  for  all  t,f  ,  (C-2) 

n 

where 

J  d(t,f )  ■  A  J'dv  exp(  1 2irvt)  D(w, f )  - 


C-l 
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These  results  agree  with  [6,  (27)4(28)].  Interpolation  rule  (C-2)  uses  the 
available  slices  of  Information  In  the  t , f  plane  of  figure  8.  A  particular 
case  of  ( C-3)  Is 

|  d(0,f )  -  (1  -  2A  |f  1 )  rect(Af)  .  (C-4) 


CAF  RECOVERY 


From  (C-l )  and  (B-7),  there  follows 

x(u,r)  -  j* df  exp(12irfr)  A.,(w,f)  0(v,f)  - 

■  X^(v,t)  ®  ^(v,t)  « 

■  xa(w,nA)  A<S(v,r  -  nA)  for  all  v,t  , 


where 


dl(v,r)  -  A  j" 


df  exp(12irf t)  D( v,f ) 


si n[,T  *(1  -  A  Ivj  )] 
it  t/A 


for  |v|  <  l 


for  all  t 


(C-5) 


0  for  \v\  >  J 

-  (1  -  A|w|)  sine  £j(l  -  A]v|)J  rect(^  Aw)  . 

Interpolation  rule  (C-5)  uses  the  available  slices  of  Information  In  the 
v.t  plane  of  figure  8.  A  special  case  of  (C-6)  Is 


(C-6) 


A^(v,0)  *  (1  -  A )v( )  rect(p  Aw)  . 


(C-7) 
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TCF  RECOVERY 


From  (C-l)  and  (B-l), 


R(t,t)  ■  £|*di»  df  exp(l2irvt  +  1 2«f t )  A^(v,f)  0(v,f)  ■ 


tT 


Mt.O  0  J’U.t) 


12 
n  1 
n+i  even 


R(^2*ia)  a2JDT< t  -  -  u)  for  an  t,t  ,  (C-8) 


whore 


a2^(t,T)  -  A2  J|*dv  df  exp(12*vt  f  12irft)  0(w,f) 


s1n2(irt/a)  -  s1n2(JJ) 

.  _  « 

4  (t2  -  4) 

A* 


sln(^)  »  s1n(**) 

}(t  *  f) 


s1n(*f)  -  »1n(g) 

St  -  b 


for  all  t,t  . 


(C-9) 


Particular  values  are 


a2£(±  j.t) 


since J)  . 


(C-10) 


It  should  be  observed  that  (C-8)  dictates  two-dimensional  Interpolation  In 
the  t,T  plane  of  figure  8.  Attempts  at  simpler  one-dimensional  Interpolation 
In  t  or  t  alone  are  bound  to  fall. 


C-3 
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PROPERTIES  OF£U,t) 


The  two-dimensional  function  £(t,T)  Is  unlike  any  proposed  previously 
for  Interpolation  of  TCP  R(t,t).  Some  of  Its  properties  are  listed  here. 

J(-t.T)  -Jftt,  -  t)  -Jtt.o  .  (c-n) 

a2 ^(iM^qA)  "  (ifll  ,  (c-i2) 

w  (m  -  d  ) 

If  m  *  q  or  m  *  -q,  then  ^'■0.  If  m  ■  q,  then  t  -  t/2,  while  If  m  -  -q, 
then  t  «  -t/2,  giving 

a2^(tqa,  2qo)  -  s1nc(2q)  .  ( C -1 3) 

Therefore JK mo, 2qa)  ■  0  for  all  m,q,  except  that  42«fr(0,0)  ■  1.  Similarly, 

& ((m  +  J)a,(2q  +  l)a)  ■  0  for  all  m.q  .  (C -14) 

If  we  define 

J(t,T)  -  ^(t,2T)  -  Sin,  CitVft),,.-  stn. .Lff-I/Al  #  (C -1 5) 

ir  (t  -  t  ) 


then 

J(T,t)  -J(1,t)  ,  (C-16) 

meaning  that  J(t,T)  Is  symmetric  about  the  45*  line  In  the  t,t  plane. 

Figure  C-l  depicts  some  sample  values  of  &  In  the  first  quadrant.  In 
particular, 
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A2  If  (mM2q  +  l)a) 


rV  -  (qt  £)2J  ’ 


d2^V(m  +  i)A,2qd)  -  -r - f 

V  2  '  ir2f(m  +  r 


C(n  +  J)2  -  q2l  ' 


(C-17) 


Interpolation  function  £(t,-r)  decays  slowest  along  the  t  -  ±t/2  lines  In 
the  t i t  plane,  and  fastest  along  the  t  ■  0  and  t  -  0  lines,  So  direct 
Interpolation  of  the  TCF  is  not  best  approximated  by  horizontal  or  vertical 
slices,  but  In  fact,  by  points  between  these  slices. 


Figure  C-1 .  Sample  Values  of  a 
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APPENDIX  0.  EVALUATION  OF  CAF 

The  original  CAF  of  waveform  s(t)  Is  given  by  (73),  In  terms  of 
spectrum  ?(f)  defined  In  (69).  As  In  (87),  we  again  presume  that  only  the 
discrete  frequency  calculations 

5(jjJ)  ■  a^>  exp(-l2irnk/N)  s(ka)  for  )n|  <  |  (0-1) 

k 

are  avallat^e.  Since  the  frequency  Increment  Is  af  ■  (NA)”1,  we 
approximate  CAF  (73)  according  to 


*(v’T)  2  na  £  #xp(12ir  NA  tJ  §(N a  f  §*<na  "  2}  " 
n 

-  f  df  txp(12«M  5(f  +  f)  5*<f  '  f)  jjj  %  ‘<f  -  jjj)  - 
J  n 

■  x(v,t)  d(t  -  n  NA)  ■ 
n 

■  5*  x(v,  t  -  n  Na)  for  all  w,t  .  (0-2) 

n 

Since  waveform  s(t)  Is  approximately  limited  to  jt|  <  T/2  (see  (84)  and 
figure  9),  then  CAF  x(»*t)  Is  approximately  limited  to  |t|  <  T,  as  may 
be  seen  by  substituting  '32)  Into  (34): 


0-1 
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x( v.t)  •  Jdt  exp(-l2m»t)  s(t  +  J)  s*(t  -  |)  .  (0-3) 

iharefore,  the  approximate  CAF  x  takes  the  appearance  shown  In  figure  D-1. 

It  Is  seen  that  overlap  Is  negligible  If  we  take 

T  <  NA  -  T,  l.e.,  N  >  .  (0-4) 

Thl  ;  requirement  on  the  F FT  size  In  (0  1)  Is  the  same  as  that  established  In 
(93)  for  the  approximate  WDF  W. 

DISCRETIZATION  IN  u  and  t 

In  order  to  utilize  available  samples  (D-1 )  In  the  evaluation  of 
app'oxlmatlon  (D-2),  we  restrict  the  evaluation  of  the  approximate  CAF  to 
frequency-shift  values 

1 exp(12’  n>  *>  ,or  a"  ’  •  <°-5> 

n 

Furthermore,  we  consider  only  the  particular  values  of  time  delay  given  by 

Kjg.  «)  •  fii  J  exp(12*rq/N)  S(=-^)  .  (0-M 

n 

since  the  right-hand  side  Is  now  an  N-poInt  FFT  for  each  m  value  of  Interest. 
N  values  of  q  are  swept  out  by  each  FFT.  Compare  (0-6)  with  (111). 


0-2 
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%(»,*) 


Figure  0-1.  Delay-Aliased  CAF 
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APPEND IX  E.  TCF  OF  IMPULSIVELY-SAMPLED  WAVEFORM 


Suppose  continuous  waveform  s(t)  Is  sampled  with  an  Infinite  Impulse 
train  (with  delay  tQ)  yielding  Impulsive  waveform 


s1(t)'  -  s(t)  i^i(t-t#-  U) 


(E-l) 


The  corresponding  TCF  Is 

^(t.T)  3  s^t  +  f)  s*(t  -  f)  - 

-  R(t,r)  “  t0  +  2  "  kA)  4(1  '  t0  “  f  “  ‘ 

k  m 

This  function  has  Impulses  In  the  t,t  plane  at 


*  -*„♦ 


2  ■  kO 


>  . 


t  -  t  -  l  «  ma 

O  L 


1  .e. ,  at 


t  ■  (k  -  m)a  , 


(E-2) 


(E-3) 


Furthermore,  the  area  of  each  of  these  Impulses  Is  1: 

JJdt  dt  4(t  -  tQ  +  j  -  ka)  4(t  -  tQ  -  g  "  mA)  ■  J*dt  4(t  +  “•<*) 

So  (E-2)  can  be  expressed  alternatively  as 


1  . 
(E-4) 


R1  (t.T)  -  R(t,t)  a22  £  mtom  H-"®  A)  4(t  “  (k  “  m)A)  •  (E“5> 

k  m 


E-1 
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Now  let 


n  »  k  +  m  ,  a  ■  k  -  m  . 
Then  nil  must  be  even,  giving 


(E-6) 


R1 (t,r)  -  R(t ,t)  A2  ^  ^  "  t0  ‘  2  A)  <(T  “  iA)  * 

n  a 
n+a  even 

-  a2  2  2.  R(t  +  f  a,  aa)  4(t  -  t  -  \  a)  «(t  -  aa  )  .  (E-7) 

n  a 
n+a  even 

This  Is  a  slight  generalization  of  (B-l),  to  allow  for  delayed  sampling. 

Thus,  the  two  approaches,  (B-1)  and  (E-l),  yield  identical  results. 

FUNDAMENTAL  TWO-DIMENSIONAL  SAMPLING  PATTERNS 

Suppose,  In  (B-l),  that  we  let 

R(t,r)  -  1  for  all  t,t  .  (E  -8) 

Then  R^t.r)  there  approaches 

ryt.t)  i  a2  ^  <(t  ••  nj)  a(x  -  aa)  .  (E-9) 

n  a 
n+J  even 

But,  at  the  same  time,  use  of  (E-8)  In  SCF  (35)  yields 

A(v,f)  ■  «(»)  4(f)  ,  (E-10) 


E-2 


k 


while  (B-5)  and  (B-8)  approach 


a,(»,n  -  2.  *(”  ■  ?)  *({  *  2i)  • 

q  m  N 
q+m  even 

Thus,  (E-9)  and  (E-11)  are  a  double  Fourier  transform  pairs 


a^(v,f) 


dt  dt  exp(-12irvt  -  I2irft)  r^t.t)  . 


(E-ll) 


(E-12) 


They  generalize  one-dimensional  result  (7)£(8)  to  two  dimensions  with 
interspersed  sampling.  The  Impulse  patterns  of  r^t.t)  and  a.,(v,f)  In 
their  respective  domains  are  displayed  In  figures  E-l  and  E-2. 
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